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ABSTRACT 

Equations  are  derived,  within  the  framework  of  geometrical 
optics,  that  relate  the  photon  flux  at  a  detector  to  the  concentration 
of  an  emitting  gas  contained  in  an  extended  source  of  rectangular  cross 
section,  for  two  optical  arrangements  that  have  been  used  experimental¬ 
ly.  Since  the  ratio  of  the  experimental  detector  signals  for  the  two 
quite  dissimilar  cases  is  in  agreement  with  that  predicted  theoretical¬ 
ly,  the  methods  are  validated. 

A  method  is  described  for  maximizing  the  signal  at  the  detector, 
using  a  given  lens,  which  may  be  applied  to  an  extended  source  of  gen¬ 
eral  cross  section. 

RESUME 

Selon  l'optique  g£om€trlque,  on  a  d£rlv£  des  Equations  qul 
relient  le  flux  de  photons  atteignant  un  d&tecteur  1  la  concentration 
d'un  gaz  €mettant  dans  une  source  Stendue  dont  la  coupe  transversale 
est  rectangulalre.  Deux  arrangements  optiques  sont  consid£r£s. 

Puisque  le  rapport  des  signaux  exp£rimentaux  pour  ces  deux  cas  tout  d 
fait  diffSrents  est  en  accord  avec  celul  prfidit  thSorlquement,  les 
m€thodes  sont  valid€es. 

Une  mSthode  est  d€crite  afin  de  maximiser  le  signal  d'un  d£tec- 
teur,  pour  une  lentille  donn&e.  Cette  mithode  peut  Stre  appliqu&e  dans 
le  cas  d 'une  source  €tendue  ayant  une  coupe  transversale  arbitraire. 
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1.0  INTRODUCTION 


During  continuing  studies  of  a  potential  gas  laser  system  it  has 
been  necessary  to  estimate  the  concentrations  of  several  energetically 
excited  molecules  in  a  vacuum  flow  system.  We  have,  for  example,  mea¬ 
sured  emissions  from  02  (alAg)  and  NF  (b*£g)  in  the  near  infrared 
(1.27  pm)  and  visible  (529  nm)  respectively,  the  spontaneous  transi¬ 
tion  rates  (kg,  s~*)  for  these  molecular  states  are  both  well  estab¬ 
lished  (1,  2),  and  so  it  is  possible,  in  principle,  to  measure  their 
concentrations  passively  in  a  gas  flow  using  calibrated  detectors  and 
narrow-band  interference  filters.  Typically,  the  detection  system  is 
arranged  to  view  the  active  gases  along  an  axis  perpendicular  to  their 
flow,  through  a  window  of  known  transmittance.  A  lens  is  frequently 
used  to  Increase  the  signal  at  the  detector.  The  field  of  view  of  the 
detector  then  covers  an  extended  source  of  variable  cross  section  con¬ 
taining  an  emitting  gas  whose  concentration  is  assumed  constant.  How¬ 
ever,  it  would  not  be  difficult  to  Include  the  case  of  a  gas  whose 
concentration  varies  in  a  known  way. 

To  calculate  an  accurate  concentration,  the  volume  of  emitting 
gas  that  can  contribute  to  the  detector  signal  and  the  fraction  of 
emitted  light  that  arrives  at  the  detector  from  any  molecule  in  this 
contributing  volume  must  be  specified  precisely. 

This  report  is  not  Intended  to  present  an  original  theoretical 
treatment  of  the  light-gathering  properties  of  optical  systems  from 
extended  sources.  This  has  been  dealt  with  elsewhere  both  in  general 
terms  (3,  4)  and  for  the  particular  case  of  optimizing  the  focal  length 
of  a  condensing  lens  to  fill  a  spectrograph  with  light  from  a  cylin¬ 
drical  extended  source  (5,  6). 
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This  document  is  concerned  with' a  somewhat  different  situation. 
It  presents  the  methods  that  have  been  used  to  calculate  accurately  the 
photon  flux  at  a  detector  under  two  optical  arrangements  employed  to 
observe  a  particular  source  cavity  that  has  a  rectangular  cross  sec¬ 
tion.  Either  a  circular  convex  lens  or  a  cylindrical  tube,  which  was 
blackened  and  baffled  to  minimize  Internal  reflection,  was  placed 
between  source  and  detector. 

However,  though  a  particular  form  of  the  extended  source  cross 
section  is  considered  explicitly,  the  methods  may  be  applied  to  ary 
cross  section.  Also,  the  treatment  of  the  lens  case  provides  a  general 

method  for  maximizing  the  signal  at  the  detector,  for  a  lens  of  a  given 

focal  length  and  diameter,  from  a  completely  general  extended  source. 

As  in  previous  work  (5,  6),  the  geometrical  optics  approximation 
is  made  (diffraction  effects  are  Ignored). 

Figure  1  defines  the  (left-handed)  Cartesian  coordinate  system 
used  in  subsequent  sections  to  describe  the  emitting  gas.  It  also 
shows  that  the  gas  is  constrained  in  the  y-direction  (±  c)  and  in  z  (0 

to  L),  but  not  in  the  x-directlon.  The  field  of  view  in  x  will  thus  be 

unconstrained. 

Since  the  mathematical  treatment  of  the  two  optical  arrangements 
(lens  and  tube)  is  quite  different,  the  theory  may  be  verified  by  com¬ 
parison  with  experiment:  the  ratio  of  experimental  detector  signals 
for  the  two  situations  (from  identical  gas  flows)  should  agree  with  the 
calculated  ratio  of  intensities  at  the  detector.  Thus  it  is  not  neces¬ 
sary  to  know  absolute  concentrations  in  order  to  verify  the  theory. 

When  this  has  been  validated,  however,  absolute  concentrations  can  then 
be  calculated  with  confidence  from  the  equations  given  here  for  either 
case.  Such  concentration  measurements  have  been  made  for  excited  oxy¬ 
gen  and  nitrogen  fluoride,  and  will  be  presented  in  separate  reports. 


UNCLASSIFIED. 

3 


Results  are  given  here  (Chapter  4)  for  oxygen  emission  intensity  ratios 
that  shew  the  agreement  between  theory  and  experiment  to  be  within  the 
measurement  errors. 

The  final  equations  for  gas  concentration  (Chapter  4)  do  not 
take  into  account  any  absorption  of  light  by  the  emitting  gas  Itself; 
l.e.  self-absorption  is  neglected.  They  are  therefore  only  valid  for 
dilute  gases,  or  transitions  for  which  the  probability  is  low,  as  in 
the  cases  of  02  and  NF  that  we  have  studied. 

This  work  was  performed  at  DREV  between  January  and  September 
1982  under  PCN  33H07,  Research  on  Chemically  Excited  Lasers. 

2 .0  LENS  CASE 


This  chapter  describes  a  method  for  calculating  the  number  of 
photons  per  second  that  arrive  at  a  circular  detector  when  a  circular 
focusing  lens  is  placed  between  it  and  the  extended  source  shown  in 
Fig.  1.  Only  thin  lens  formulae  will  be  used. 

The  emitting  volume  will  be  treated  as  a  continuum  of  contribu¬ 
ting  object  planes  (xy)  to  be  mapped  onto  the  detector  plane.  The 
fraction  of  each  mapped  surface  Intercepted  by  the  detector  will  be 
derived.  This,  together  with  a  function  specifying  the  solid  angle 
subtended  to  the  lens  by  each  point  In  the  object  plane,  will  be  inte¬ 
grated  through  z-space  to  yield  a  total  photon  flux  at  the  detector  for 
a  given  gas  concentration  and  emission  rate . 

Originally,  the  cylindrical  symmetry  imposed  by  the  coaxial 
circular  lens/detector  arrangement  will  be  exploited,  the  actual  cavity 
dimensions  being  treated  as  constraints  upon  this. 
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2.1  The  General  Contributing  Object-Plane  Radius 

Figure  2  shows  a  cross  section  through  the  xz  plane  of  Fig.  1. 
A  lens  of  focal  length  f  and  radius  R  is  placed  at  a  distance  s  from 
the  window  (at  z  •  0)  of  the  extended  source.  The  detector  (radius  r) 
is  a  distance  v  from  the  lens.  The  bar  above  the  v,  as  with  other 
symbols  used  here,  indicates  a  fixed  parameter:  v  is  a  particular 
value  in  the  complete  Image  space  v. 

There  exists  a  u  (in  the  object  space  u)  for  a  fixed  (v,  s) 
configuration,  such  that  an  object  of  radius  h  is  mapped  exactly  into 
an  image  of  radius  r  at  the  detector  position,  u  is  given  by 

u  ”  fv/  (v  -  f)  [  l] 

provided  v  >  f,  and  u  lies  between  s  and  s  +  L. 

The  object  radius,  ¥  is 


h  ■  u  .  r/v  [2] 

The  general  object  position  is 

u(z)  -  Z  +  S  [3] 

with  u  being  a  particular  value  z  +  s. 

At  any  other  position  on  the  z-axls  (z  ■  0,  L),  the  maximum 
object  radius  that  can  contribute  to  the  intensity  at  the  detector  is 
given  by  hf  or  h-  (see  Fig.  2).  This  is  because  only  rays  passing 
through  both  h  and  the  lens  can  arrive  at  the  detector.  This  is  not  to 
say  that  all  rays  leaving  any  h(u)  arrive  at  the  detector,  as  we  shall 
see. 
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So  far,  no  way  has  been  specified  for  choosing  the  lens  and 
detector  positions  v  and  s.  Conventionally,  the  object  position  u  is 
required  to  be  at  the  centre  of  the  cavity  (z  -  L/2),  and  the  object 
radius  h  is  chosen  to  be  the  cavity  radius.  Equations  1  and  2  are  then 
implicitly  simultaneous  in  two  unknowns  (v,  s).  However,  for  the  case 
considered  here  in  which  the  source  is  unrestricted  in  the  x-direction, 
h  is  not  clearly  defined.  For  the  moment,  [l]  and  [2]  simply  define  71 
and  h  for  a  given  choice  of  v  and  s. 

It  is  necessary  to  have  an  expression  for  the  general  contribu¬ 
ting  object  radius  h(u): 

a)  in  u  <  u  ,  h  =  h-  in  Fig.  2  and  the  geometry  gives: 

(R  -  h-)/u  -  (R  -  h)/u  [4] 

using  [2]  yields: 

h-  ■  r.u/v  +  R(u  -  u)/u  [5] 

b)  in  u  >  u  ,  h  =  h+ 

(R  +  h+)/u  -  (R  +  h)/u  [6] 

and  similarly, 

h+  ■  r.u/v  +  R(u  -  u)/u.  [7] 

Thus  in  general,  for  L  +  a  >  u  >  s: 

h(u)  ■  r.u/v  +  r|u  -  u|/u  [8] 

Note  that  h(u)  reduces  to  h  (eq.  2)  as  required. 
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h(u)  defines  a  continuum  of  objects  that  are  radii  of  cross 
sections  In  the  xy  plane.  Each  of  these  objects  will  be  focused,  pro¬ 
vided  s  >  f,  at  a  different  position  on  the  image  side  of  the  lens. 
Therefore  it  is  necessary  to  know  the  radius  of  the  circle  at  any  point 
in  the  Image  space  through  which  all  rays  from  any  object  circle  of 
radius  h(u)  pass,  so  that  the  fraction  of  this  light  intercepted  by  the 
detector  can  be  specified. 

2.2  Light  Cone  Radius  at  any  Position  on  the  Image  Side  of  the  Lens 

Consider  Fig.  3.  An  object  of  radius  h  at  u  Is  focused  at  v 
with  Image  radius  la.  The  radius  of  the  cross  section  of  the  cone 
through  which  all  light  from  the  object  passes  Is: 

1-  at  v-  <  v,  and 
1+  at  v+  >  v. 


Now, 


tan9  (R  -  i„)/v  -  (R  -  i-)/v-  [9] 


so  that 


1-  ■  i0  v-/v  +  R(v  -  v-)/v  [lO] 

Now  consider  the  detector  to  be  placed  at  v-;  i.e.  v-  =  v. 


With 


i0  ■  vh/u 


[11] 
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u  v 


FIGURE  3  -  Light  cone  radii  on  the  image  side  of  the  lens 

Eq.  10  becomes 

i-  -  hv/u  +  R(v  -  v)/v  [ 12 ] 

for  v  <  v. 

Similarly  for  v+  >  v 

tan*  -  (R  +  i+)/v+  -  (R  +  i„)/v  [  1 3] 

which  leads  to 

i+  -  hv/u  +  R(v  -  v)/v  [  1 4] 

for  v  >  v. 

Thus,  at  a  fixed  position  v  (detector  position),  the  radius  of 
the  circle  of  light  from  an  object  of  radius  h  may  be  written  in  a 
general  form,  using  [l2]  and  [l4]  as  a  function  of  object  position  u 
only: 


i(u)  ■  h(u)v/u  +  R  v  -  v  /v 


[15] 
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where  h(u)  is  given  by  [8]  for  all  u,  and 


v(u)  ■  fu/(u-f) 


Using  [8]  and  [l6],  i(u)  can  be  written  explicitly 


i(u)  -  r  +  2Rv 


I  -  I 

U  ~  u 


which  may  also  be  written  in  the  useful  form: 


i(u)  ■  hv/u  +  Rv 


1  _  1 
u  u 


2 .3  Object  and  Image  Space  Cross  Sections 


[16] 

as  a  function  of  u: 

[17] 


[18] 


Equation  8  for  h(u)  defines  a  continuum  of  circular  emitting 
cross  sections,  and  [l8]  for  i(u)  gives  the  circles  into  which  these 
are  mapped  at  the  detector  position  v  The  y-direction  constraint 
imposed  by  the  cavity  dimension  ±  c,  (see  Fig.  1)  can  now  be  included. 


The  shaded  region  of  Fig.  4  shows  the  constrained  emitting  cross 
section,  whose  area  is: 

A„(z)  ■  2 [ h2  sin“l  (c/h)  +  c(h2  -  c2)1/2]  [ 19] 


for  h  >  c. 

Figure  5  shows  the  cross  section  A^(z)  into  which  A»(z)  is 
mapped.  Its  value  is  given  by: 


Ai(z)  “  2t1x  8ln_l  (1y/ix)  +  S  (1x  "  ty1'2]  t2°] 
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for  >  1^,  which  results  from  h  >  c.  lx  may  be  calculated  from  [l7] 
since  It  Is  generated  by  h: 


(*) 


r  +  2Rv 


1 

Q 


1 

u 


[21] 


whereas  iy,  since  it  results  from  c,  must  be  calculated  from  [  18] 


iy(z)  -  cv/u  +  Rv 


[22] 


The  fraction  of  A^(z)  that  Is  intercepted  by  the  detector  can 
now  be  specified. 


2.4  Fraction  Intercepted  by  Circular  Detector 

A1(z)  may.  In  general,  cover  the  circular  detector  (of  radius  r) 
In  four  ways: 


(a)  1  >  r  and  1  >  r  (Fig.  6(a)). 

y  x 

In  this  case  the  fraction  of  A^(z)  that  Is  received  by  the 
detector  Is: 

Fi(z)  -  itr2/A1(z)  [23] 


A^(z)  being  given  by  eqs.  20  to  22. 
(b)  ly  >  r  ;  ix  <  r. 
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FIGURE  6  -  Interception  of  A^(s)  by  a  circular  detector 
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This  case  does  not  occur  because  It  implies  1  >  1  ,  which  can- 

y  x 

not  occur  if  h  is  chosen  to  be  everywhere  >  c  (by  adjustment  of  the 
distances  v  and  s). 

(c)  i  >  r  ;  i  <  r  (Fig.  6  (c)). 

^  y 

The  shaded  area  of  Fig.  6(c)  is 

Ad(z)  -  2[r2  sin"1  (iy/r)  +  iy  (r2  -  i2)l/2],  [24] 

and 

Ft(a)  -  Ab(z)/A1(z)  [25] 

(d)  i  <  r;  i  <  r  (Fig.  6(d)). 

*  y 

In  this  case  the  entire  area  A^z)  is  intercepted  by  the  detec¬ 
tor,  and 

Fj(z)  -  1  [26] 


To  summarize: 

(1)  at  all  points  in  z,  we  have  defined  an  emitting  area  A«(z) 

(eq.  19)  that  can  contribute  to  light  intensity  at  the 
detector; 

(2)  through  eqs.  23  to  26,  we  have  obtained  the  fraction  of  this 
area  that  is  mapped,  by  the  lens,  onto  the  detector;  and 

(3)  it  remains  to  define  the  fraction  of  light  emitted  by  any 
point  in  A«(z)  that  arrives  at  the  lens  (essentially  the 


T 


UNCLASSIFIED 

13 


solid  angle  subtended  to  the  lens),  and  hence  by  (1)  and  (2) 
above,  the  fraction  of  emitted  light  that  arrives  at  the 
detector. 


2 .5  Fraction  of  Emitted  Light  Arriving  at  the  Lens 


This  calculation  is  complicated  by  the  fact  that  in  the  y- 
direction  light  emission  is  restricted  by  an  opaque  roof,  and  that  to 
support  a  window,  the  cavity  has  a  wall  of  thickness  T  where  there  is 
essentially  no  emitting  gas. 
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window  lens 


FIGURE  8  -  Area  of  lens  Illuminated:  AL(z) 

The  case  Y  >  R  occurs  when  the  radiating  point  Is  close  enough 
to  the  window  to  Illuminate  the  entire  lens  area. 

Thus  the  solid  angle  subtended  to  the  lens  (A^/u2)  falls  off 
more  rapidly  than  simply  1/u2,  since  the  cavity  roof  restricts  the  area 
of  the  lens  that  can  Intercept  light;  l.e.  Aj^  ♦  0  as  z  Increases. 

The  fraction  of  light  emitted  by  a  point  in  A,(z)  that 
arrives  at  the  lens  is  given  by: 
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Fl<«)  "  AI(z)/Airu2  [31] 

Equation  31  is  exact  only  for  points  on  the  z-axis,  but  for  small  off- 
axis  variations  it  is  a  very  good  approximation. 

2.6  Relation  Between  Emitting  Gas  Concentration  and  Photon  Flux  at  a 
Detector 


If  the  radiative  rate  constant  for  spontaneous  emission  from  the 
gas  is  ks  (s_1),  and  the  gas  concentration  is  a  (molecules  cm“3), 
then  an  infinitesimal  volume  element 

6V  -  A„(z)  6z  (cm3)  [32] 


emits 


6P  -  k8aA0(z)6z  (photons  s-1)  [33] 


in  all  directions. 

Of  these,  F^(z).dP(z)  arrive  at  the  lens,  and  Fi(z).  F^z). 
6P(z)  arrive  at  the  detector. 

Integrating  over  the  extended  volume  (0  to  L  in  z-space)  gives 
the  total  number  of  photons  per  second  arriving  at  the  detector, 

PSL>  ignoring  for  the  moment  the  transmittance  of  the  window,  the 
lens  and  a  narrow-band  interference  filter.  Thus: 

PSL  ”  o  Fi<*>V‘)ks°A*(*)d*  (Photons  s_l)  [34] 


or,  from  [31 ] 
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vsLm-hr0  Vz)AL(z)Ae(z),u'2*dz 


[35] 


The  subscripts  on  P  indicate  spontaneous  and  lens  respectively. 

It  is  convenient  to  define  a  lens  integral: 

1^  =  /  F^zjA^CzjAoCz)  .u“2  .dz  [36] 

so  that  [35]  can  be  written  in  a  form  that  will  be  useful  in  future 
comparisons  with  results  obtained  for  the  viewing  tube  case.  Equa¬ 
tion  35  becomes: 


V  ■  v-v4'  l371 

2.7  Calculation  of  the  Lens  Integral:  A  Numerical  Example 

A  FORTRAN  program  has  been  written  to  calculate  the  individual 
terms  in  the  Integrand  of  eq.  36  pointwise  in  z-space,  and  then  per¬ 
form  the  integral  numerically  using  Simpson's  rule.  The  precision  of 
the  result  can  always  be  verified  by  increasing  the  number  of  integra¬ 
tion  points.  The  following  dimensions  are  required  as  input  parame¬ 
ters: 

r  -  active  detector  radius; 

R  -  lens  radius; 

f  -  lens  focal  length; 

c  -  cavity  y-dimension  (Fig.  1); 

T  -  cavity  wall  thickness  (Fig.  7); 

L  -  cavity  z-dlmension  (Fig.  1); 
s  -  lens-to-source  distance;  and 
v  -  lens-to-detector  distance. 
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Column  one  of  Table  I  lists  a  particular  set  of  values  {qj}  for  these 
parameters,  which  correspond  to  an  experiment  using  a  gas  flow  con¬ 
taining  excited  oxygen  ((^(a^A)).  The  emission  from  this  (at  1.27  pm) 
was  recorded  with  a  Judson  J-16  germanium  detector  through  a  narrow- 
band  interference  filter. 

The  value  of  the  integral  for  this  input  set,  IjJq^},  was 
found  to  converge  to  five  significant  figures  (0.42401  cm3)  using  201 
integration  points. 

Table  I  also  gives  an  error  analysis.  The  experimental  uncer¬ 
tainties  {Sq*}  in  the  measured  dimensions  {qjJ  are  listed  in 
column  two,  and  the  percentage  uncertainties  in  column  three.  The 
total  measurement  uncertainty  is  thus  about  8Z.  Column  four  shows  the 
integral  values  obtained  using  a  particular  parameter  at  its  maximum 
allowed  value  (qj  +  6qj)  with  all  other  parameters  held  at  the 
values  given  in  column  one.  Thus,  column  five  shows  the  maximum  per¬ 
centage  variation  in  1^  that  can  be  introduced  by  the  uncertainty  in 
each  of  the  input  parameters.  The  total  uncertainty  in  1^  is  then 
about  9Z  due  tc  the  cumulative  measurement  error.  Hence: 

I  -  0.4240  ±  9Z  (cm3)  [38] 

Li 

2.8  Choice  of  v  and  s:  Maximization  of  IL 


In  the  absence  of  any  other  criteria,  the  lens  and  detector 
positions  of  Table  I  (s  and  v)  were  chosen  such  that  an  object  at  the 
centre  of  the  cavity  (z  ■  L/2),  with  a  radius  of  about  one  and  a  half 
times  the  cavity  height  (l.e.  h  ■  1.5  c),  would  form  an  image  at  the 
detector  with  a  radius  equal  to  that  of  the  detector  (r).  This  ensures 
h(z)  >  c  for  all  z  in  0  to  L,  as  required  by  eq.  19. 
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TABLE  I 

Input  and  error  analysis  In  the  calculation  of  1^ 


Parameter  value 

qi 

Uncertainty 

6qi 

X  6qt 

IL(qj+iqJ) 

X6I 

L 

r 

0.412 

0.001 

0.2 

0.42585 

mm 

R 

2.23 

0.01 

0.4 

0.42586 

f 

6.35 

0.05 

0.8 

0.41897 

c 

0.54 

0.02 

3.7 

0.44325 

T 

1.21 

0.01 

0.8 

0.42364 

0.09 

L 

12.05 

0.05 

0.4 

0.42486 

0.20 

S 

12.9 

0.10 

0.8 

0.42675 

0.65 

V 

9.5 

0.10 

1.0 

0.41941 

1.08 

(q  In  centimeters;  I  -  cm3) 
X  L 


However,  it  is  possible  to  optimize  s  and  v  by  numerically  maxi 
mixing  the  lens  Integral  1^  of  eq.  36,  and  hence  by  [37]  ensuring 
that  the  photon  flux  at  the  detector  is  a  maximum. 


The  numerical  procedure  is  as  follows:  a  particular  s-value, 
Sj,  is  defined  (Sj  >  f);  1^  is  calculated  for  a  series  of  v  values 
converging  to  a  maximum  1^  (e.g.  using  the  Newton-Raphson  method)  and 
hence  also  to  an  optimum  Vj.  A  pair  (vj»  8j)  can  thus  be  defined  for 
any  s  >  f  that  maximizes  1^,  and  by  scanning  through  s  an  absolute 
maximum  1^  and  its  associated  optimum  pair  (v,  s)  can  be  obtained. 


Figures  9  and  10  show  the  results  of  such  a  calculation  for  the 
experimental  setup  described  in  Section  2.7.  The  absolute  maximum  of 
1^  occurs  at: 
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8  *  22.1  ;  v  "  8.14  (cm)  [39] 

Equations  1  and  2  can  be  used  to  calculated  the  size  (h)  and 
position  (u)  of  an  object  that  is  imaged  exactly  on  the  detector,  for 
the  optimized  values  of  [39].  This  yields: 

u  -  28.9 

z  “  u  -  s  “  6.8  [40] 

and  h  -  1.46. 

The  nominal  object  position  (z)  is  thus  displaced  from  the 
centre  of  the  cavity  (at  z  ■  6)  and  the  object  radius  (h)  is  almost 
three  times  the  cavity  height. 

These  results  are,  of  course,  peculiar  to  the  extended  source 
considered  here,  which  is  restricted  in  two  dimensions  and  unbounded  in 
the  third,  and  also  to  the  lens  and  detector  dimensions.  However,  the 
methods  would  be  applicable  to  any  extended  source  for  which  it  is 
important  to  maximize  a  detector  signal. 

3.0  TUBE  CASE 


The  equations  derived  in  this  chapter  specify  the  number  of 
photons  per  second  arriving  at  a  circular  detector  through  a  cylindri¬ 
cal  tube  placed  between  the  extended  source  of  Fig.  1  and  the  detector. 
The  tube  was  blackened  and  baffled  to  minimize  internal  reflection. 

The  purpose  of  this  viewing  tube  is  to  define  clearly  the  contributing 
volume  of  emitting  gas. 

Figure  11  shows  the  xz-plane  of  Fig.  1,  indicating  a  tube  of 
length  K  and  radius  R  (exaggerated  for  clarity),  and  the  volume  of  gas 
that  can  contribute  to  light  Intensity  :t  the  detector. 
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Any  point  in  the  volume  bounded  by  h(z)  (for  z  ■  0  to  L)  in 
Fig.  11  may  illuminate  the  entire  detector  surface.  However,  in  the 
region  h(z)  to  H(z)  (shaded  in  Fig.  11),  radiating  molecules  can  only 
illuminate  a  fraction  of  the  detector  area.  This  fraction  goes  from  1 
to  0  as  x  goes  from  h  to  H. 

A  first  approximation,  considering  only  the  volume  bounded  by 
h(z),  will  be  presented,  and  subsequently  the  exact  equations  will  be 
derived  for  the  entire  volume  bounded  by  H(z).  It  will  be  shown  that 
the  contribution  from  the  shaded  volume  of  Fig.  11  is  not  at  all  negli¬ 
gible;  in  fact  it  is  of  the  same  order  of  magnitude  as  the  first 
approximation  (for  the  tube  length  considered  here).  Clearly,  it  will 
go  to  zero  as  the  tube  length  increases,  but  unfortunately  so  will  the 
detector  signal.  In  practice  the  tube  length  is  therefore  limited  and 
the  exact  treatment  must  be  used. 

3.1  First  Approximation:  Detector  Wholly  Illuminated 
From  the  geometry  of  Fig.  11: 

h(z)  -  u(z).  +  r  [4l] 


where,  as  in  chapter  2, 


u(z)  ■  z  +  s  [42] 

An  emitting  area,  Aj(z),  in  which  any  point  can  fully  illuminate 
the  detector,  can  thus  be  defined.  Figure  12  shows  this  cross  section, 
bounded  by  the  cavity  height  (c)  and  the  variable  h(z). 


It  is  helpful  in  understanding  the  subsequent  treatment  of  the 
shaded  region  of  Fig.  11,  to  derive  the  area  A^(z)  in  cylindrical  polar 
coordinates: 
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h  ir/2  p,  ,  .  , 

Ax(z)  -  4{/x  /  pdpd#  +  /  /*  pdpd*}  [43] 

*1 

where  #j  and  Plt  shown  in  Fig.  12,  are  given  by 

^(z)  -  sin-1  (c/h(z) )  [44] 

and 

Pl(®)  ■  c  cosec  *  [45] 

Substitution  of  [44]  and  [45]  into  [43]  leads  to  the  familiar 
form  (cf.  eq,  19)  for  the  area  of  a  section  of  a  circle: 


Aj(z)  •  2(h2  sin-*  (c/h)  +  c(h2  -  c2)^2)  ;  (h  >  C)  [46a] 

which  may  also  be  written: 

A^(z)  -  2(h2#1  +  c2cot  *j)  [46b] 

(if  c  >  h,  Aj (z)  -  ifh2) 

An  infinitesimal  volume  element 

5V  ■  pfip6#6z  [47] 

centred  at  a  position  u(z),  containing  a  molecules  cm-3,  would  emit 
k^afiV  photons  per  second  into  a  sphere  of  surface  area  4i»u2.  If  all 

points  in  6V  can  fully  illuminate  a  detector  of  area 


a  ■  nr 


2 


[48] 


then  the  number  of  photons  per  second  that  arrive  at  the  detector  from 
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^I’TlV67  f49] 

4x  u*  8 

Integrating  over  the  volune  bounded  by  h(z)  and  the  cavity 
length,  and  using  [43]  and  [46]: 

P,  -  — —  /  dz  [50] 

1  4ir  °  u2 

which  is  valid  provided  h(z)  is  small  relative  to  u(z).  Using  [48]: 

k  or2 

P1--V-1!  [51] 


which  defines  It  (cf.  [ 50 ] ) . 

Px  is  thus  a  first  approximation  to  the  number  of  photons  s_1 
arriving  at  the  detector  from  the  cavity,  which  ignores  contributions 
from  the  region  h  to  H  of  Fig.  11.  In  that  region,  an  emitting 
molecule  illuminates  only  a  fraction  of  the  detector  area. 

3.2  Inclusion  of  Boundary  Region 

Suppose  that  an  infinitesimal  volume  element  (eq.  47)  contains 
molecules  that  irradiate  a  fractional  area  of  the  detector: 

a(p,z)  -  F(p,z).wr2  [52] 

F(p,z)  is  the  fraction  of  the  total  area  (nr2)  that  is  irradi¬ 
ated  by  a  molecule  at  (p,  4,  z),  thus: 

0  <  F(p,z)  <  1  [53] 

This  fraction  is  independent  of  4  (cylindrical  symmetry),  but  is 
dependent  on  p,  and  through  the  limits  of  p  (h  and  H),  on  z. 
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Equation  49  must  now  be  replaced  by: 

5pst  ■  ^  •  Vsv  t54) 

Using  [52]  and  integrating  over  the  volume  bounded  by  the  cavity  dimen¬ 
sions  and  H(z)  now, 


ks“r2  v  F(p  z) 

P  =  * —  ;  dv' 

ST  4  „2 


[55] 


and  by  analogy  with  [43]: 


k“r2  L 


H 


c  ^  f  ^9  tt/2  p  1  .  r  . 

rST  "  “V“  f°  [4/o  d*  l  F<P.Z>  OdP  +  {  d*  /1F(P,r)Pdp]  [56] 
u  *2 


The  subscripts  on  P  indicate  spontaneous  and  tube  respectively. 

The  limits  (#2  an<^  H)  of  the  integrals  in  [56]  now  include  the 
region  h  to  H.  From  the  geometry  of  Fig.  11: 


H(z)  -  u(z)  .  +--r^  ~  r 


[57] 


and  by  analogy  with  [44]  (cf.  Fig.  12  with  h  replaced  by  H,  and  by 
♦2): 


*2(z)  -  sin"1  (c/H(z) )  [58] 

Now,  if  a  weighted  area  A(z)  is  defined: 

A(z)  =  4{/2  ;  F(p,z)pdpd*  +  }/2  l1  F(p,z)pdpd*}  [59] 

*2 

then  [56]  can  be  written  in  a  form  similar  to  [50]  and  [ 51 ] : 


i 
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PST  “ 


k  or2 
8 


A  (z) 


dz  = 


k  ar2 
8 


[60] 


which  defines  1^,  a  tube  integral  (cf.  the  lens  Integral  case, 
eq.  3  7). 

A(z)  is  a  weighted  area  in  that  points  in  an  area  A(z)  (Fig.  12 
with  h  replaced  by  H,  and  *x  by  #2)  are  weighted  according  to  the  frac¬ 
tion  of  the  detector  they  are  capable  of  illuminating. 


Equation  59  can  be  separated  into  the  previously  defined  area 
A}(z)  (eq.  43),  and  a  weighted  area  Ag (z)  associated  with  the  region  h 
to  H  in  p-space.  In  0  <  p  <  h,  F  ■  1,  thus  [59]  becomes: 


_  .  *.>H  *,h 

A(z)  -  4{/  zf  F  pdpd*  +  /z/  pdpd* 

#  h  #  • 

*/2  Pi  4 1  pi  h  .  r  _ 

+  /  Pdpd*  +  f1  fl  Fpdpd*  +  f1  f  pdpd*]  [ 6l] 

*1  *2  h  *2  ° 


which,  by  comparison  with  [43]  can  be  seen  to  yield: 

A(z)  -  4[*2(z)  /  F(p ,z)pdp  +  f1  j1  F(p,z)pdpd*}  +  Aj(z)  [62] 


or 


A(z)  =  A2(z)  +  Ax(z) 


[63] 


defining  A2(z). 

Thus  it  is  possible  to  compare  the  effect  of  including  the 
shaded  area  of  Fig.  1  with  the  first  approximation  given  by  eqs.  50, 
51.  Equation  60  with  [63]  gives: 
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P 


ST 


k  or2 

-T-l*. 


[64] 


l.e. 


k  ar2 

PST*-24—  t1!  +  ^  t65l 


where 


I2  =  \  •  dz  [661 

u 

It  only  remains  to  define  F(p,z),  the  fraction  of  the  detector 
area  Illuminated  by  a  molecule  at  p  between  h(z)  and  H(z).  This  is 
algebraically  complicated,  and  is  treated  In  detail  In  Appendix  A. 
Conceptually,  It  is  only  necessary  to  know  that  an  explicit  algebraic 
form  exists  for  F(p,z)  in  terms  of  tube  and  detector  dimensions  and  the 
distance  s,  which  satisfies  [53],  so  that  the  integrals  In  [62]  can  be 
calculated  numerically. 

It  Is  convenient  for  such  a  calculation,  to  write  F  in  the 

form: 


F(p,z)  -  1  -  G(p, z)  ;  (h  <  p  <  H;  0  <  z  <  L)  [67] 

where  G  varies  from  0  to  1  as  p  goes  from  h  to  H. 

Equation  62  may  then  be  expressed,  after  some  calculus  and 
algebra,  as: 

A(z)  -  4{«2(z)[*(H2  -  h2)  -  /  G(p,z)pdp]  -  %h2(*l  -  *2) 

h 

-  fl  ?l  G(p,z)pdpd*  +|  [(H2  -  c2)*  -  (h2  -  c2)^]]+A1(z)  [68] 

♦2  h  L 

\ 

j 


j 
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Defining  the  two  Integrals  in  [68]  to  be  Ip  and  Ip#,  and  then  rearrang¬ 
ing: 

A(z)  -  4{%[*2H2  -  *^2]  -  *2Ip  -  Ip#  +|  (cot*2  -  cot*!)} 

+Aj(z)  [69] 

Aj(z)  can  thus  be  eliminated  by  using  [46b] : 

A(z)  -  4{*2(^H2  -  Ip)  +  %c2  cot*2  -  Ip#}  [70] 

Equation  70  is  a  rather  simple  form  for  directly  calculating  the 
weighted  area  A(z)  needed  for  the  tube  integral  1^  of  eq.  60.  The 

integrals  in  p-  and  4-space  are: 

r  (*)  s  /  G(p,z)pdp  [71] 

“  h 

I  (z)  5  f1  d*  /  G(p,z)pdp  72 

p*  *2  h 

If  c  >  h,  eq.  70  is  still  valid,  but  with  *!  replaced  by  w/2  in 
[72]  for  Ipi. 

3.3  Calculation  of  the  Tube  Integral:  A  Numerical  Example 

A  FORTRAN  program  calculates  both  the  first  approximation  inte¬ 
gral  Ix  (eq.  51)  and  the  accurate  form  1^  (eq.  60)  using  [46]  and  [ 70] 

for  the  numerators  in  the  integrands.  The  z-space  Integrations  are 
performed  numerically  using  Simpson's  rule.  For  1^,  at  each  point  in 
z-space  a  p-space  Integral  [7l]  is  performed  for  Ip(z),  and  a  double 
Integral  is  performed  for  Ip^(z)  (eq.  72):  at  each  point  In  *-apace  a 

p-space  integral  mist  be  calculated  since  the  upper  limit  of  the  p~ 
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Integral  is  9-dependent.  All  the  functions  vary  smoothly,  and  a  rapid 
convergence  to  five  significant  figures  is  obtained  with  only  25  points 
for  each  numerical  integration. 

The  required  input  dimensions  are: 

L  -  cavity  z-dimension  (Figs.  1  and  11) 
c  -  cavity  y-dimension  (Fig.  1) 

8  -  detector-to-source  distance  (Fig.  11) 

R  -  tube  radius  (Fig.  11) 
r  -  detector  radius  (Fig.  11) 

K  -  tube  length  (Fig.  11) 

Table  II  lists  a  set  of  experimental  values  for  these  parame¬ 
ters,  together  with  an  error  analysis  as  in  Table  I,  corresponding  to 
an  experiment  using  a  gas  flow  of  excited  oxygen  as  described  in  Sec¬ 
tion  2.  7. 

The  total  measurement  error  for  this  input  is  about  8%,  which  is 
shown  in  Table  II  to  give  a  comparable  total  uncertainty  in  the  tube 
Integral  Ij.  Thus: 


IT  -  0.4135  E-01  ±  8%  (cm.)  [73] 

The  first  approximation  Ij  was  0.2160E-01,  which  is  only  52%  of 
Ij.  Thus  the  contribution  to  the  photon  flux  at  the  detector  from 
the  boundary  region  (h  to  H,  the  shaded  region  of  Fig.  11)  is  of  a 
similar  magnitude  to  that  from  the  region  in  which  all  molecules  can 
fully  illuminate  the  detector.  In  general,  the  full  integral  oust  be 
calculated. 


I 
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TABLE  II 


Input  and  error  analysis  in  the  calculation  of  If 


Parameter  value 

«i 

Uncertainty 

6qt 

X  &qt 

lOx 

Lj^q^+dqj) 

%6lT 

L 

12.05 

0.05 

0.4 

0.41484 

0.33 

c 

0.54 

0.02 

3.7 

0.42648 

3.14 

s 

18.2 

0.20 

1.1 

0.41028 

0.77 

R 

0.49 

0.01 

2.0 

0.42371 

2.47 

r 

0.412 

0.001 

0.2 

0.41341 

0.02 

K 

11.9 

0.1 

0.8 

0.40973 

0.91 

(all  values  in  centimeters) 


4.0  LENS/ TUBE  RATIO:  THEORY  AND  EXPERIMENT 


For  the  cavity  of  Fig.  1  and  a  lens-detector  setup  as  in  Fig.  2, 
eq.  37  gave  the  number  of  photons  per  second  arriving  at  the  detector 
This  must  be  modified  to  include  transmission  factors  for  the  window 
(wtf),  lens  (ltf)  and  narrow-band  Interference  filter  (ftf),  thus: 

PgL  -  (wtf).(ltf).(ftf).ksaIL/4Tr.  [74] 


For  example,  a  quartz  lens  has  92Z  transmission  at  1.27  um,  in 
which  case  ltf  “  0.92. 


For  the  same  cavity,  but  with  a  tube-detector  setup  (Fig.  11), 
the  photon  flux  at  the  detector  was  given  by  [ 60 ] .  Again  modified  to 
include  transmission  factors  (no  lens)  this  becomes: 

PST  -  (wtf).(ftf)ksaIT.r2/4  [75] 

Thus  for  a  given  gas  emitting  at  a  rate  kg,  with  a  fixed  con¬ 
centration  a,  viewed  with  the  same  detector  through  the  same  filter, 
the  ratio  of  the  lens  and  tube  arrangements  is: 
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(ltf) 

ir2 


[76] 


Experimentally,  the  detector  signals  are  measured  as  voltages 
(Vi,  and  V<r  for  lens  and  tube  cases  respectively),  which  are  related 
to  photons  s-1  by  a  calibration  constant  T,  thus: 

p„„  -  rvT  ;  P„  -  H?  [  77] 

SL  L  ST  T  1  J 


Hence  an  experimental  value  for  the  ratio  of  lens  to  tube  inte¬ 
grals  can  be  obtained  from  [77]  and  [76]: 


-  »r. 

I  (ltf) 

expt. 


[78] 


For  the  excited  oxygen  emission  experimental  setups  considered 
in  Sections  2.7  and  3.3,  theoretical  values  for  1^  and  Ij  were 
obtained  (eqs.  38  and  73),  giving  the  theoretical  ratio: 


0.4240 

IT  “  0.04135 
theory 


10.3 


[79] 


Table  III  shows  the  results  from  a  series  of  experiments  in 
which  different  total  oxygen  pressures  were  used  in  the  cavity,  i.e. 
different  values  of  a,  the  emitting  02  concentration.  The  correspond¬ 
ing  detector  signals,  measured  using  the  aforementioned  tube  and  lens 
arrangements,  and  their  ratios  are  given.  The  ratio  obtained  at  any 
fixed  pressure  should  be  Independent  of  the  pressure.  The  table  shows 
that  this  is  the  case  since  no  ratio  deviates  by  more  than  2%  from  the 
average  (18.83). 
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TABLE  III 


Detector  readings  for  several  gas  pressures 


Pressure 

(torr) 

VL 

(yv) 

vT 

(uV) 

VL/VT 

1 

185 

9.7 

19.1 

2 

32  7 

17.7 

18.5 

3 

443 

23.7 

18.7 

4 

535 

- —  — 

28.0 

19.1 

The  experimental  value  for  the  ratio  I^rlf  then  given  by 
[78],  with  r  -  0.412  and  ltf  ■  0.92, 


XL  _  it  (0 .412)2 

I_  “  02 

T 

expt. 


(18.83) 


10.9 


[80] 


The  difference  between  the  theoretical  and  experimental  ratios 
is  less  than  6%,  which  is  well  within  the  uncertainty  introduced  into 
the  theoretical  ratio  (about  172!)  by  uncertainties  in  the  measured 
dimensions. 


The  absolute  concentration  a  of  an  emitting  gas  in  the  cavity  of 
Fig.  1  is  related  to  the  measured  voltage  from  a  calibrated  detector: 

a)  in  the  lens  case  using  eqs.  74  and  77, 

o  -  4*rvL/(wtf)(ltf)(ftf)kgIL  ;  [ 8 1 ] 

b)  in  the  tube  case  using  [75]  and  [77], 

a  -  4rvT/(wtf)(ftf)kflIIx2  .  [82] 
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The  integrals  Ij,  and  If  are  calculated  from  physical  dimen¬ 
sions  only  (Chapters  2  and  3),  and  T  is  the  detector  calibration  con¬ 
stant.  The  integrals  may  be  computed  to  a  high  precision,  but  the 
values  are  only  as  accurate  as  the  total  dimension  measurement  error, 
which  is  typically  less  than  10%.  The  accuracy  of  the  calculated  con¬ 
centration  is  also  dependent  upon  the  accuracy  of  the  known  emission 
rate  (k8)  which  is  about  10%  in  the  excited  oxygen  case,  and  the 
detector  calibration  (<  10%).  The  transmission  factors  and  voltage 
measurements  can  be  quite  accurately  measured  (about  1%) . 

An  overall  concentration  measurement  that  is  accurate  to  within 
about  30%  can  therefore  be  expected. 

Absolute  concentrations  of  a  particular  excited  species,  02(8*4) 
in  the  gas  flow  from  a  chemical  generator,  have  been  estimated  by 
observing  two  entirely  different  emissions  arising  from  this  molecule. 
The  first  of  these  emissions  (at  1268  no),  the  direct  spontaneous 
transition  to  the  ground  state,  was  observed  with  a  germanium 
detector/filter  combination  using  both  tube  and  lens  arrangements,  and 
gave  rise  to  the  voltages  listed  in  Table  III.  From  these,  partial 
pressures  of  02(a*A)  were  calculated  in  the  range  from  0.22  to 
0.63  torr  as  the  total  02  pressure  varied  between  1  and  4  torr.  Obser¬ 
vation  of  the  second  emission  (the  red  "dimole"  emission  at  634  nm) 
using  a  silicon  detector  and  appropriate  filter,  led  to  estimated 
02(a*A)  partial  pressures  between  0.27  and  0.69  torr  over  the  same 
total  pressure  range. 

Agreement  between  the  estimates  from  the  two  emissions  is  within 
the  anticipated  30%  accuracy,  and  the  average  of  the  calculated  per¬ 
centages  of  02(a1A)  in  the  flow  (varying  between  25  and  17%  over  the  1 
to  4  torr  total  02  pressure  range)  is  of  the  expected  magnitude  for 
chemically  generated  singlet  delta  oxygen  (Ref.  7)  in  the  flow  system 
employed  in  these  experiments. 
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5.0  CONCLUSIONS 

The  concentration  of  a  photon-emitting  gas  contained  in  an 
extended  source  can  be  estimated  using  a  calibrated  detector,  narrow¬ 
band  Interference  filter,  and  either  a  lens  or  a  viewing  tube  to  define 
the  contributing  gas  volume. 

The  close  agreement  between  theory  and  experiment  for  the  ratio 
of  intensities  at  a  detector  for  the  two  cases  gives  us  confidence  in 
the  methods  and  formulae  developed  in  this  report.  The  final  two  equa¬ 
tions  could,  of  course,  both  be  in  error  by  the  same  multiplicative 
constant.  Though  specifically  applied  to  a  source  of  rectangular  cross 
section,  the  methods  can  be  modified  for  any  extended  source  without 
difficulty,  and  the  lens  treatment  provides  a  method  for  optimising  the 
light  intensity  at  a  detector  from  such  a  general  source. 
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As  the  point  moves  off-axis  and  into  the  region  where  h  <  p  <  H, 
this  projected  circle  essentially  moves  across  the  detector.  There 
will  be  some  distortion  from  a  circle,  and  Rp  will  increase  slightly 
as  p  becomes  large.  However,  for  small  off-axis  displacements  the 
projection  can  be  considered  to  be  a  displaced  circle  of  radius  Rp. 

It  would  not  be  difficult  to  write  the  exact  eliptical  equation  for  the 
projected  curve,  but  the  deviation  from  this  is  negligible  here 
because:  a)  we  are  concerned  only  with  the  area  of  intersection  of  the 
detector  circle  and  the  projection,  the  deviation  from  the  circular 
curvature  of  which  is  minimal  in  the  small  arc  covering  the  detector; 
and  b)  at  the  extreme  value  p  ■  H,  where  the  deviation  from  a  circle  is 
the  greatest,  this  is  also  where  the  contribution  to  the  integral  IT 
Is  the  smallest  because  the  fractional  area  illuminated  is  close  to 
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z 

►  P 


FIGURE  A-I  -  Illumination  of  a  circular  detector  through  a  cylindrical 
tube 


N>|  =»  Nj|  =< 
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FIGURE  A-2  -  Area  of  intersection  of  two  circles 


Now,  the  area  of  intersection  of  two  circles  (Fig.  A-2)  of  radii 
rj  <  r2  is  given  by: 

r $  -  (a  -  £)[r$  -  (o  -  O2]^  ~  r|  sin-1  (°  ^1) 

-  S(r?  -  K2)h  -  r?  sin"l  (|-)  [A-3] 

where  o  is  the  distance  between  the  centres  of  the  two  circles,  and 

5  -  [a2  -  (rg  -  rf)]/2o  [a-4] 

In  the  case  considered  here, 
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r2  =  Rp  ,  rj  =  r  [a-5] 

and  o  may  be  calculated  from  the  geometry  of  Fig.  A-l: 

a(z, p)  -  Kp/(u-K)  [A-6] 

Hence  F(p,z),  the  fractional  area  of  the  detector  illuminated 
is: 

F(p,z)  -  A  /irr2  -  /r)2  -  { (a-?) [r2-( o-£)2] *  +  R28in-1(^.)}/irr2 

Ip  p  P  Kp 

+  *  -  €  (r2  -  52)Vwr2  -  [sin-^i)]/*  [A-7] 

Equation  A-7  may  be  rearranged  to  give  a  form  (cf.  eq.  67): 

F(p,z)  -  1  -  G(p,z)  [A-8] 

where 

G(p , z)  -  [T(p,z)/ir  -  (R2  -  r2)/2]/r2  [a-9] 

and 

T(P,z)  -  (o  -  0[R2  -  (o  -  O2]*  +  R2  sin"1  (HlL) 

P  P  Rp 

+  C(r2  -  C2)^  +  r2  sin-1  (i)  [a-10] 

Thus,  equations  A-8  to  A-10,  together  with  [a-6]  for  o,  [a-4] 
for  S,  and  [ A— 2 ]  for  Rp,  define  the  fractional  area  required. 
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